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A  TECHNIQUE  FOR  THE  SYNTHESIS  OF  LINEAR,  NONSTATIONARY 

FEEDBACK  SYSTEMS 
PART  I.  THE  APPROXIMATION  PROBLEM 


A.  R.  Stubberud 

Summary:  In  this  paper  and  in  a  companion  paper  a  technique  for  the  synthesis  of  linear, 
nonstationary  feedback  systems  is  presented.  This  paper  deals  with  the  approximation 
problem,  i.e. ,  a  method  is  developed  whereby  a  realizable  weighting  function  which 
produces  a  desired  output  when  subjected  to  a  polynomial  input  can  be  generated. 

Certain  conditions  which  must  be  fulfilled  by  the  outputs  are  also  developed. 

INTRODUCTION 


Literally,  the  phrase  "synthesis  of  a  system"  implies  that  a  designer  proceeds 
logically  in  a  step-by-step  manner  from  a  set  of  system  specifications  to  a  system 
which  will  meet  these  specifications.  In  this  sense  many  of  the  design  techniques  for 
control  systems  are  not  true  synthesis  techniques.  As  examples,  in  the  field  of  linear, 
stationary  control  systems,  the  root  locus  and  frequency  domain  methods  are  not  true 
synthesis  techniques.  On  the  other  hand,  the  design  method  proposed  by  Guillemin  1, 
in  which  the  closed  loop  transfer  function  is  first  determined  from  the  specifications 
and  the  compensation  is  determined  subsequently  is  a  synthesis  technique. 

Synthesis  techniques  in  the  areas  of  nonlinear  control  and  linear,  nonstationary 
control  are  few  and,  in  general,  quite  specialized.  No  general  synthesis  techniques 
for  feedback  systems  in  either  of  these  areas  have  been  developed.  The  reason  for 
this  is  quite  simple.  No  general  analytical  techniques  comparable  to  the  techniques 
available  for  stationary  systems  are  available  for  nonlinear  and  linear,  nonstationary 
systems.  The  analytical  techniques  which  are  available  are  usually  quite  difficult, 
quite  specialized  or  both  and  therefore  do  not  lend  themselves  readily  to  application  to 
the  synthesis  problem. 

Engineers  have  put  more  effort  into  developing  general  approaches  to  the  analy¬ 
sis  of  linear,  nonstationary  systems  than  they  have  into  nonlinear  systems.  The  two 
general  approaches  which  have  received  the  most  attention  are  the  weighting  function 
approach  2,  3,  4  and  the  time  varying  transfer  function  approach.  5,  6,  7  Each  of  these 
is  a  logical  extension  of  stationary  system  analysis.  Another  approach  to  analysis 
which  has  been  discussed  to  a  lesser  degree  is  the  differential  equation  approach. 

8,  9,  10,  11 

Synthesis  techniques  for  nonstationary  systems  have  tended  less  toward  certain 
general  approaches  than  have  the  analytical  techniques.  One  approach  which  has  been 
proposed  synthesizes  a  given  nonstationary  transfer  function,  or  a  given  weighting 
function,  as  an  open  loop  system  consisting  of  parallel  elements  of  lower  order.  12, 13 


Another  synthesis  scheme  allows  a  given  nonstationary  weighting  function  to  be  synthe¬ 
sized  as  an  open  loop  system  composed  of  analog  computer  elements.  2,  4,  14, 15  Pro¬ 
cedures  which  allow  a  nonstationary  system  to  be  synthesized  in  the  form  of  a  feed¬ 
back  system  are  relatively  few.  Malchikov,  16,  proposed  one  procedure  by  means  of 
which  a  given  weighting  function  can  be  synthesized  as  a  feedback  system.  The  pro¬ 
cedure  is  applicable  to  the  case  where  the  system  is  required  to  contain  a  "fixed" 
plant.  This  method  has  the  disadvantage  that  certain  integral  equations  must  be  solved 
at  specific  points  in  the  development. 

One  possible  way  of  categorizing  the  techniques  of  analysis  and  synthesis  of 
linear,  nonstationary  systems  is  to  classify  them  according  to  the  mathematical  repre¬ 
sentations  which  are  used  to  describe  the  elements  of  the  system.  In  this  sense  the 
three  categories  would  be:  (1)  nonstationary  transfer  functions;  (2)  weighting  functions; 
and  (3)  differential  equations.  These  categories  cannot  be  made  too  rigid  since  in  any 
one  problem  it  might  be  desirable  to  use  more  than  one  representation.  Each  of  the 
three  categories  has  its  own  methods  for  combining  elements  of  a  system  and  each  has 
certain  advantages  over  the  others.  As  was  indicated  in  the  preceding  paragraph, 
synthesis  techniques  which  fall  into  the  first  two  categories  have  been  developed;  how¬ 
ever  none  have  been  developed  which  fall  into  the  third.  The  differential  equation 
approach,  however,  has  certain  combinatorial  techniques  which  are  advantageous  in 
the  synthesis  problem  and  the  purpose  of  this  paper  is  to  develop  a  synthesis  technique 
which  exploits  these  combinatorial  techniques. 

SCOPE  OF  THIS  PAPER  AND  COMPANION  PAPER  18 


Guillemin,  1,  in  1947,  proposed  that  the  synthesis  of  linear,  stationary  feedback 
systems  should  be  accomplished  in  three  steps  as  follows: 

(a)  The  closed-loop  transfer  function  is  determined  from  the  specifications. 

(b)  The  corresponding  open-loop  transfer  function  is  found. 

(c)  The  appropriate  compensation  networks  are  synthesized. 

In  this  paper  this  synthesis  technique  is  generalized  to  include  linear,  non¬ 
stationary  systems  which  may,  or  may  not,  be  required  to  contain  a  "fixed"  plant. 
This  technique  makes  use  of  the  combinatorial  techniques  developed  in  Reference  8. 
The  steps  to  be  followed  in  the  nonstationary  case  are: 

(a)  The  closed  loop  weighting  function  is  determined  from  the  specifications. 

(It  is  desirable  to  specify  the  system  by  its  weighting  function  since  the 
system  output  for  any  input  can  then  be  determined  by  convolution.  ) 

(b)  From  the  closed  loop  weighting  function  determine  the  corresponding 
closed  loop  differential  equation. 

(c)  Determine  the  open  loop  differential  equation  from  the  closed  loop  differ¬ 
ential  equation. 
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(d)  Determine  the  differential  equations  for  the  appropriate  compensation  net¬ 
works  and  synthesize  these  networks  by  appropriate  analog  computer  systems. 

Implementation  of  this  procedure  requires  techniques  for:  (1)  determining  a 
weighting  function  which  meets  the  given  specifications;  and  (2)  manipulating  this 
weighting  function  in  such  a  way  that  the  appropriate  compensation  networks  can  be  de¬ 
veloped.  A  method  is  developed  by  means  of  which  a  weighting  function  can  be  de¬ 
termined  which  meets  the  specifications  that:  (1)  the  input  to  the  system  is  a  polynomial 
in  time;  and  (2)  the  corresponding  output  approximates  a  given  time  function.  In  ad¬ 
dition  the  application  of  the  combinatorial  manipulations  developed  in  Reference  8  to 
the  synthesis  problem  are  described  in  detail  in  a  related  companion  paper.  18  To 
clarify  these  procedures  some  simple  examples  are  worked  in  this  and  the  companion 
paper. 18 


THE  APPROXIMATION  PROBLEM 

The  first  step  of  the  synthesis  problem  is  to  determine  a  system  weighting 
function  which  will  perform  according  to  the  given  specifications.  This  step  will  be 
called  the  approximation  problem.  Actually  for  the  technique  which  is  developed  in 
this  report  it  would  be  adequate  to  determine  a  system  differential  equation;  however  it 
is  more  desirable  to  know  the  system  weighting  function  since  the  system  output  can 
then  be  determined  for  any  input  by  application  of  the  convolution  integral. 


There  are  many  methods  by  which  a  system  weighting  function  might  be  obtained; 
for  instance,  it  might  be  given  directly  in  the  system  specifications.  In  any  event,  once 
it  has  been  established  (it  must  have  certain  special  properties  defined  below),  the  pro¬ 
cedure  presented  in  the  companion  paper  18  can  be  followed  to  complete  the  synthesis 
problem.  In  the  rest  of  this  section  a  specialized  method  for  determining  weighting 
functions  is  presented.  The  method  is  restricted  to  systems  which  receive  inputs 
which  are  expressed  as  polynomials  in  time  and  whose  outputs  can  be  approximated  as 
a  separable  function  (see  below). 


This  approximation  method  can  be  formulated  as  follows.  In  Figure  1  let 


Xlt) 


w 
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LINEAR  SYSTEM 
FIGURE  1 


x(t)  be  a  polynomial  in  time  given  by 

N 

x(t)  =  x(t  -t)  =  /  a  (t-r)*^  t^T 

=  0  t<T  (1) 

where  the  a  are  constants,  t  is  time,  and  t  is  the  time  of  application  of  x(t)  to  the 
linear  system  W.  y(t)  is  the  output  of  the  linear  system  W  which  is  produced  as  a  re¬ 
sult  of  the  application  of  x(t).  It  is  assumed  that  an  analytic  expression  for  y(t)  can  be 
determined  from  the  specifications.  W  is  the  unknown  linear  system  whose  weighting 
function  W(t,T)  is  to  be  specified  in  a  form  which  can  be  synthesized  as  a  feedback 
system.  From  the  convolution  integral,  y(t)  can  be  expressed  as 
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(2) 


w(t,  e)x{e)dd 

-00 


Substituting  Equation  (1)  into  Equation  (2),  y(t)  becomes 

^  t  r  i"' 

y(t)  =  y^(t,T)  =  ^  ®n  y  W(t,e)|^  e-rjde;  t 


n«o 


=  0 


t<T 


(3) 


The  problem  now  reduces  to  that  of  solving  the  integral  Equation  (3)  for  W(t,T).  Fortu¬ 
nately  this  equation  reduces  readily  to  a  linear,  stationary  differential  equation  which  is 
easily  solved.  Prior  to  showing  this,  however,  it  is  necessary  to  restrict  the  class  of 
W(t, t)  which  are  allowed  in  the  synthesis  technique.  Similar  restrictions  are  placed  on 
y  (t, t)  later.  It  is  assumed  in  all  the  succeeding  works  that  W(t,T)  corresponds  to  the 
solution  of  the  following  linear  differential  equation 


(4) 


where  6(t-T)  is  a  delta  fvmction  occurring  at  time  t,  and  c^  (t)  and  b^  (t)  are  continuous 

functions  of  time.  The  properties  of  a  W(t,  t)  satisfying  an  equation  with  the  form 
of  Equation  (4)  have  been  discussed  in  detail  previously  2,  4,  8, 10, 11, 12, 13, 14  and 
this  discussion  will  not  be  repeated  here.  Suffice  it  to  say  that  W(t,T)  can  be  written  in 
the  form 


W(t,T)  =  W,(t,T)  +  b  (t)  6(t-T)  t  St 

1  m 

=  0  t  <T 


(5) 


where  Wj(t,T)  has  the  form 

m 


Wj(t,T)  = 

J=1 


(t) 


=  0 


t  St 


t<T 


(6) 


in  which  the  q's  must  have  2n  continuous  derivatives  and  the  ^'s  must  have  n  continuous 
derivatives.  8 

Substituting  Equation  (5)  into  Equation  (3),  the  relationship 
N  .  ^  N 


^  f  W^(t,T)[  9-T]"d9  +  b^(t)  i  a^(t-r)" 
n=o  T  n=o 


(7) 
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results.  The  solution  of  this  equation  will  be  broken  into  two  parts:  (1)  determination 
of  b  (t),  and  (2)  determination  of  W,(t, t). 


First,  a  technique  for  determining  b  (t)  will  be  presented.  Assume  that  a, 

m  k 

(0  2  k  S  N)  is  that  one  of  the  non-zero  a^'s  with  the  lowest  valued  subscript,  i.  e. , 


N  ^  N 

y  (t,T)  =  /  a  f  W,  (t,e)(e-T)“d0  +b_(t)  ^  a  (t-T)" 
X  n  1  m  n 


N 

I 


(8) 


If  now  the  kth  partial  derivative  of  with  respect  to  t  is  formed,  the  relationship 


d^y^(tjT) 


(t,  0)(-l)Nn)(n-l). .  .  (n+k+l)(e-T)”'‘^d0  + 


+  b  (t)a  (-l)^(n)(n-l). . .  (n-k+l)(t-T)"  ^ 
m  K 


(9) 


9 

results.  Now  if  the  limit  of  - r -  as  T-*t  is  formed,  the  equation 

3t 


9 


8t 


=  b  (t)  a,  (-1)  (k!) 
m  k  '  ' 


(10) 


T-*t 


is  formed.  Solving  Equation  (10)  for  b  (t),  the  desired  relationship  for  b  (t)  is 

m  m 

obtained,  i.e.. 


(-1) 


"  (k!)a 


T-^t 


(11) 


Once  b  (t)  is  known  the  expression 
m 


N 

b  (t)  /  a  (t-r)^ 
m  ^  n 
n=o 


(12) 


can  then  be  formed.  Now  a  new  function  y^'(t,T)  is  formed  from  Equation  (8)  as 
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(13) 


N  N  ^ 

y  ’(t.T)  =  y  (t.T)  -b^(t)  y  a^(t-T)''  =y  a  f  W, (t.  0)(0-T)''d0 
X  X  m  Lj  n  Lj  nj  i 

n=o  n®o  T 

The  second  problem  now  is  to  determine  Wj(t,T)  from  Equation  (13).  First,  deriva¬ 
tives  of  y^'(t,T)  with  respect  to  t  are  formed.  The  first  derivative  is: 


9yjj'(t.T) 

9t 


N  ^ 

=  ^  an  f  Wj(t,e)(-l)n(e-T)'''^ 

n=o 


de  =  a  W,(t,T) 
o  1 


(14) 


The  kth  derivative  (k<N  +  1)  is: 


8^y^'(t,  t)  ^ 


dr 


t 

"  Z  ^n  j*  Wj(t,e)(-l)V)(n-l). 


n»o 

k-1 


-1 

3=0 


!) 


8^’^’^Wj(t,T) 


8t 


k-l-j 


(n-k+l)(e-T)"'^d0 


k  =  1,2,.  ..,  N 


(15) 


j,N+l  .  N 

9  y..(t.T)  V' 


8^'^W,(t,T) 


8t 


N+1 


V'  o  w,v 

j=o  •’  8t 


(16) 


Equation  (16)  is  seen  to  be  a  linear,  stationary  differential  equation  whose  solution  is 
Wj(t,T).  To  solve  this  differential  equation  uniquely  N  initial  conditions  must  be  ob¬ 
tained.  Examination  of  Equations  (15)  shows  that  in  these  equations  if  r-^t. 


8  y^'(t,T) 
8t‘^ 


k-1 


=  ^  (-l)^‘*‘\j!)a. 

T-t 


8*^‘^'^Wj(t,T) 


8t 


k-l-j 


(17) 


T-*t 


or 

/ 

8^‘V(t,T) 

-■J  ' 

.  O 

8^y^'(t,T) 

ar"-* 

dr^ 

T-^t 

\ 

k-2 


T-*-t 


j+1  (t,r) 

'T-*t 

k  =  1,2,  ... ,  N 


(18) 
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Equations  (18)  represent  a  set  of  recursive  relationships  which  provide  the  N  necessary 
"initial  conditions".  Actually  since  the  derivatives  of  Wj(t,T)  are  with  respect  to  t. 

Equations  (18)  represent  final  conditions.  This  inconvenience  is  easily  remedied  by 
making  the  substitution 

T  =  t  -  z  (19) 

into  Equations  (16)  and  (18).  The  resulting  equations  can  then  be  easily  solved  and  by 
making  the  inverse  substitution 

z  =  t  -  T  (20) 

into  the  solution  obtained,  Wj(t,T)  is  recovered. 

As  a  special  case  of  this  development  Equations  (11)  and  (16)  are  examined  for  the 
degenerate  case  of  x(t),  i.e., 

x(t-T)  =  aj^(t-T)^  (21) 

Equation  (11)  becomes 


b 

m 


(t) 


(-U 


gy  ^(t,T) 
N®  X 


8t 


N 


T  =  t 


(N!)  a 


N 


Equation  (16)  becomes: 


(22) 


j,N+l  . 


3t 


N+1 


(-l)^'^^(N!)a„  W,  (t,T) 
N  1 


(23) 


or 


(-1) 


N-fl 

N+1  9”^  y^'(t.T) 


Wj(t,T) 


8t 


N+1 


(N !  )a 


N 


(24) 


Equations  (23)  and  (24)  then  should  be  used  in  the  event  that  the  input  is  a  polynomial 
with  but  one  non-zero  coefficient.  This  case  obviously  includes  step  function  inputs. 


It  is  obvious  that  in  order  for  W(t,T)  to  have  the  form  indicated  in  Equations  (5)  and 
(6)  certain  restrictions  must  be  placed  on  the  form  of  Referring  to  Equation  (1) 

let  X  (t-r)  *  a  (t-r)",  then  x(t-T)  can  be  rewritten 
n  n 
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(25) 


N 


N 


X(t-T)  =  ^  a  (t-T)"  =2 

n*o  n*o 


Equation  (3)  can  then  be  rewritten  in  the  form 

N  ^ 

y^(t,T)=^  f  W(t.0)  x^(0)  de 
*  n=o  -00 


(26) 


However,  since  W(t,T)  satisfies  Equation  (4),  yjj(t.T)  and  x(t-T)  are  related 
through  the  equation: 


m  -i  .  m  .  m  N 

y  b.(t,%i>  =  y  Yb. 

i=o  J=o  *^“0 


(t) 


a^X  (t-T) 

n 


at 


j 


(27) 


N  m 


PO  1=0 


1-  a^x  (t-r) 

y  b  (t)  ^ 


Examination  of  the  input  x(t-T)  term  by  term  reveals  that  each  x^(t-T)  =  a^(t-T) 
can  be  formed  by  integrating  the  delta  function  (a^)(n!)  6(t-T)a  total  of  n+1  times,  i.e 


n 


X  (t-r)  =  a  (t-r)^  =  C- -^-C a  n!  6(t-T)  dt 
n  n  J  J  n 


(28) 


Using  relationships  (26),  (27)  and  (28)  a  block  diagram  which  corresponds  to  Equation 
(27)  can  be  drawn  as  shown  in  Figure  2. 


In  Figure  2,  y  is  that  part  of  the  output  y(t)  due  to  the  x  (t-r)  component  of 
V  .  n  n 

x(t-T),  i.e., 

N 


y  (t,T)  =  /  y„(t,T) 

X  n 


(29) 


n=o 


The  notation  1  /p*'  within  a  box  in  Figure  2  indicates  that  the  input  to  the  box  is  integrated 
n  times .  The  system  in  Figure  2  apparently  represents  a  linear  system  with  a  delta 
function  as  an  input  and  y  (t,T)  as  an  output.  Figure  3(a)  represents  the  nth  branch  of 
the  system  in  Figure  2. 
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1 


I 


BLOCK  DIAGRAM  EQUIVALENT  OF  EQUATION  (27) 

FIGURE  2 


•  BRANCH  OF  FIGURE  2 

FIGURE  3 

» 
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yx(t,T) 


1 


Each  branch  then  is  seen  to  consist  of  a  cascade  combination  of  two  linear  differential 
equations.  This  combination  can  be  replaced  by  a  single  differential  equation  (see  the 
companion  paper  18)  as  indicated  in  Figure  3(b).  Figure  2  is  then  equivalent  to  N+1 
parallel  differential  equations  and  it  can  be  replaced  by  a  single  equivalent  linear  differ¬ 
ential  equation  (see  companion  paper  18).  Since  the  input  to  this  system  is  a  delta 
function  and  the  output  is  y  (t.r),  then  yx(t»T)  is  obviously  the  weighting  fimction  of  the 
equivalent  linear  differential  equation.  The  first  requirement  on  y  (t,T)  then  is  that  it 
must  satisfy  all  of  the  requirements  of  a  weighting  function,  i.  e. ,  u  must  be  a  separa¬ 
ble  function  of  t  and  t  with  the  form 

M 

y  (t,T)=y  Qf  (T)u  (t)  t^ 

jfe'o  J  J 

=  0  t<T  (30) 

where  the  a.  (t)  are  linearly  independent  and  the  u.(t)  are  linearly  independent.  Obvi- 
J  3 

ously  yx'('t«T)  (see  Equation  13)  will  have  the  same  form.  The  possibility  of  y  (t,  t) 
containing  singularity  functions  (i.e. ,  delta  functions  or  derivatives  of  delta  functions) 
is  dismissed  since  in  practice  other  than  continuous  outputs  will  not  be  encountered. 

In  addition  to  the  requirement  that  y  (t,  t)  have  the  form  indicated  in  Equation  (30), 
examination  of  Equations  (11),  (16),  and  (18),  which  are  rewritten  here  for  convenience, 
reveal  some  additional  requirements  on  y^(t,T). 


b^(t) 


(-1)^  8t^ 


(k!)a, 


(11) 


9t  ]=o 


1)'’  (j!)a 

J 


a^'^Wj(t,T) 


9t 


N-j 


3^'Vl(t,T) 


9t 


k-1 


T-»t 


a 


k-1 


9^  '(t,T) 

X _ 

9t^ 


T-*t 


j=i 


T-^t 


k  =  1,2,  ..  .,N 


(16) 


(18) 
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The  solution  for  Equation  (16)  has  the  form 

N 


W 


Y  -■I'jd-T)  pT  ct) 

j(t,T)>2v  G(a-T)(-1)  — rw-’*" 

i=l  t  da 


(31) 


where  the  c.  are  functions  of  the  "initial  conditions"  in  Equation  (18);  the  6’s  are  roots 
of  the  polynomial 


N 

^  (j ! )  a  .7^'^  =  0 


and  G(x-6)  is  the  solution  of  the  differential  equation 


3=0 

From  Equation  (30) 


M 


y  '(t,  a)  =  /  a  (a)  u  (t) 
*  3  3 


and 


i=l 


.N+l  V 

9  yjj'it.T) 


N+l 


M 

=  1 


( a)  u^(t) 


9o‘-  -  i^l 

If  G(a-T)  is  written  in  its  general  exponential  form  as 

N 


G( a-r) 


I 


T  e 
r 


-y^ict-T) 


(32) 


(33) 


(34) 


(35) 


(36) 


then  Wj(t,T)  can  be  written  in  the  form 


i=l  i=l  t  j=l 


(CT)u.(t)dCT  (37) 
3 


Rearranging  Equation  (37) 

N 


^  y.rr  -7.t  ^  1 

Wj(t,T)  =  2  ®  '  +(-1)^'^^  T.  2  u.(t)F..(t)J 

3=1 


1=1 

M 


M  N 

■f  u.(t)[(-l)”2  T.e’^'V  (T)] 

3=1  i=l 


(38) 
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where 


(39) 


(A  -y^a 

F^j(t)  •  J  e  Qfj  (of)  dff 

Equation  (38)  indicates  that  W^(t,T)  is  a  separable  function  of  t  and  r  with  a  maximum 
of  M+N  terms;  therefore  it  can  be  written 

M+N 

Wi(t,T)  g.(T)u.(t)  (40) 

j=l 

where  the  g's  and  u's  are  defined  according  to  Equation  (38).  In  order  for  Wj(t,T)  to  be 

a  weighting  function  of  the  type  considered  in  Equation  (6)  it  is  sufficient  that  the  g's 
have  M+N  continuous  derivatives  and  the  u's  have  2(M+N)  continuous  derivatives. 

Applying  these  differentiability  requirements  to  Equation  (38)  it  can  be  seen  that  the  a's 
should  have  2M+3N  continuous  derivatives  and  the  u's  should  have  2M+2N  continuous 
derivatives. 


Example 

As  an  example  of  the  technique  just  developed,  consider  the  following  problem. 

The  input  to  a  linear  system  has  the  form: 

x(t-T)  =  2u(t-T)  +  (t-r)  -  2(t-T)^  (41) 

From  the  system  specifications,  the  output  of  the  system  is  required  to  have  the  form: 


.t^ 

■  5  ^  10  1 

'[2T+ 

1  2 
2*^ 

2  3./ 

+  gT  +1 

+  t 

_-3T^+iT®  + 

ir+T+z] 

(42) 


From  the  given  input  x(t-T)  and  the  desired  output  the  weighting  function  W(t,  t),  which 
represents  the  desired  linear  system,  is  to  be  determined.  It  was  shown  that  W(t,  r)  can, 
in  general,  be  divided  into  two  parts,  i.e., 

W(t,T)  =  W,(t,T)  +  b  (t)  6(t-T)  (43) 

1  m 


The  output  can  then  be  represented  by 


N  ^  N 

y  (t,T)  =y  a  f  W,(t,e)(e-T)''d0  +b  (t)  y  a  (t-T)*' 
•'x  /U  n  J  1  m  ^  n 

n=o  T  n-o 


(44) 
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where  N  *  2,  a  =  2,  a  «  1,  and  a_  =  -2. 

Ox  A 

The  first  step  in  determining  W(t,T)  is  to  find  b  (t).  Taking  the  limit  of 
+  m 

y^(t,  t)  as  t-»T  ,  the  relationship  in  Equation  (45)  is  obtained. 


lim  y  (t,T)  »  2t  =  b  (t)  a  • 
*  m 


t-^T 


(45) 


therefore 


b  (t)  =  ^  =  t 
m  a 

o 


(46) 


The  second  term  of  in  Equation  (44)  is  then 


1 

N 

*  tr2-T-2T^J  +  t^ 

n*o 

(47) 

1 

1 

1 

Inserting  (47)  into  (44)  and  solving  for 

I 

1 

N 

y.,'(t.T)=  y  (t.T)  -  b  (t)  )  a  (t-T)" 

X  X  m  n 

n*o 

(48) 

* 

Equation  (49)  is  obtained: 

<  • 

I/*  X  ^  *4r5  .  10  1^.3 

y,'(t,T)  =  -5t  +t  L-+-Tj+t 

[s-fr  -T^] 

x2ro  12  231  J 

+  t  [-2T+2T  +3T  J  +  t| 

2^1  3^1  4] 

.  +6^  ^6^  J 

(49) 

Next  the  derivatives  of  y^'(t,  t)  with  respect  to  t  are  formed 


a^x'  .  10 

ar  “  3 


t^  +  t^[-  I  -  2t]  +  t^  [-2  +  T  +  2t^  ]  +  t[-2T  +  I 


2  2  3 

T  +  gT  _ 


(50) 


aV  ' 

•'x 


-2t^  +  t^  1  +  4t  J  +  t  -2  +  T  +  2t^  j 


(51) 


aV  '  2 

- =  4t^  +  t  +  4tT 

dr 


(52) 
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1 


Referring  to  Equations  (16)  and  (18)  in  the  previous  section  it  can  be  seen  that 


9  y  '(t,T)  2  aw,(t,T) 

- 3 -  *  4t  +  t  +  4tT  *  -2a2Wj(t.T)  +  a^  - 


-  a 


8t 


Wj(t,t  ) 


1 

a 


a^^'(t.T) 


8t 


T  -►t 


9  Wj(t.T) 
°  8t^ 


2t 


9W^(t,T) 

8t 


2 

1 

S  M  MM  / 

9  y  ' 

-aW(t,t')i 

s 

^  a 

T  -►t  o 

3t 

k 

T-t 

/ 

>  =t 


The  solution  of  the  set  of  Equations  (53),  (54)  and  (55)  then  is  the  desired  weighting 
function.  To  simplify  the  algebra  involved  and  to  make  the  solution  of  this  set  of 
equations  easier,  the  following  substitutions  are  made 

Z  =  t  -  T 

k2  =  2t 
ki  =  4t 


(53) 


(54) 


(55) 


f(z)  =  Wj(t,  t-z) 


(56) 


Then  since  a  =  2,  a  =  1  and  a-  *  -2,  Equations  (53),  (54)  and  (55)  can  be  rewritten  as: 

OX  A 


d  f  ^  1  df  .  . 

— 7  +  o  —  2f  =  k_z  -  k, 
,  2  2  dz  2  1 

dz 


(53') 


f(0)  =  2t 
f'(0)  =  -t 

The  complementary  solution  of  Equation  (53')  has  the  form 

-Tl*  -7,* 

f  (x)  =  Ae  +  Be 
c 


(54') 

(55') 

(57) 


where  A  and  B  are  fvuictions  of  the  initial  conditions  and  7^  and  are  roots  of  the 
polynomial 

2  1  o  « 

7  +  j  7-2  =  0 


(58) 
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The  particular  solution  has  the  form 

k_ 

2  2 


*  2  8  2 

The  entire  solution  then  has  the  form 

-7,x  -7„x  k  k-  k„ 

f(x)  =  Ae  ^  +  Be  ^  ^  ^  ^  x 

Applying  the  initial  conditions  of  (54')  and  (55'),  the  constants  A  and  B  can  be  deter¬ 
mined  to  be 

k_  k.  n  r  1  k„ 

f(0)  + 


(59) 


(60) 


^2  1  r  1  ^2 


(61) 


1  -2 


B  = 


^2  r  ^2 

f.(0)  +  ^+7^Lf(0)+/-^J 


^1  -^2 


(62) 


Substituting  the  values  in  (54'),  (55')  and  (56)  it  can  be  shown  that 
A  =  0 
B  =  0 

f(x)  therefore  equals  the  particular  solution  of  the  differential  equation,  i.  e. , 

^2  ^2 

f(x)  =  ^  -  8“  ■  y  * 


(63) 


(64) 


Substituting  the  values  of  Equations  (56)  into  (64)  it  is  seen  that 

Wj(t,  t)  =  t^  +  tr  (65) 

The  overall  weighting  function  defined  by  Equation  (43)  is  then 

W(t,T)  =  t^  +  tT  +  t6(t-T)  (66) 

and  the  problem  is  solved.  Convolution  of  (66)  and  (41)  will  indeed  produce 
as  given  in  Equation  (42). 

APPROXIMATION  OF  SEPARABLE  FUNCTIONS 


From  the  discussion  in  the  previous  section  it  is  obvious  that  for  a  large  class  of 
synthesis  problems  the  response  function  must  have  the  form  of  a  separable  function 
which  satisfies  suitable  differentiability  conditions.  In  addition,  the  class  of  weighting 
f;mctions  considered  in  this  paper  are  separable  functions.  It  is  then  necessary  to  the 
synthesis  technique  to  have  a  method  of  approximating  a  fimction  of  two  variables  as  a 
separable  function  of  the  two  variables.  Some  work  has  been  done  on  this  problem;  for 
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instance,  Cruz,  17,  proposed  an  impulse-train  approximation  of  weighting  functions;  and 
Cruz  and  Van  Valkenberg,  13,  proposed  a  method  of  approximating  weighting  functions  as 
separable  functions  which  approximate  the  weighting  functions  in  some  sense.  The  rest 
of  this  section  describes  a  simple  method  of  approximating  a  function  of  two  variables  as 
a  separable  function. 

As  a  tool  to  help  explain  this  approximation  method,  a  final  value  controller  will 
be  considered.  Suppose  that  a  system  is  to  be  designed  according  to  the  following  speci¬ 
fications.  The  system  will  receive  a  step  input  at  some  time  (which  is  unknown)  between 
t  s  0  and  t  =  T  and  it  is  required  that  the  system  settle  to  within  n%  of  its  final  steady 

state  value  (unity)  by  time  T  where  T  >  T.  On  the  other  hand,  it  is  desirable  to  keep 

s  s 

the  system  as  sluggish  as  possible  at  all  times.  Thus  if  the  system  is  designed  as  a 
nonstationary  linear  system,  the  first  step  of  the  design  is  to  determine  a  step  response 
y(t,T)  of  this  system  which  satisfies  the  specifications.  From  this  y(t, t),  the  appropriate 
weighting  f;mction  can  be  determined  and  the  design  completed. 

Figure  4  represents  a  two-dimensional  view  of  y(t,T)  in  which  the  t  axis  is  per¬ 
pendicular  to  the  plane  of  the  paper.  The  planes  y(t,  0)  and  y(t,  T)  are  the  two  planes 
which  bound  the  region  of  interest  of  y(t,T),  i.e. ,  if  the  function  y(t,  t)  is  suitable  in  this 


TWO  DIMENSIONAL  VIEW  OF  y(t,  T) 


FIGURE  4 


region  the  specifications  can  be  met. 
form: 


y(t,  0)  =  1  -  e 
y(t,  T)  =  1  -  e 


-oit 


-a2 


(t-T) 


Suppose  that  y(t,  0)  and  y(t,  T)  are  chosen  with  the 

(67) 


where  a  and  a„ 

X  M 


are  chosen  so  that 


y(T  ,  0)  =  y(t  ,  T)  =  1 
s  s 


n 

100 


(68) 


It  is  seen  then  that  y(t, t)  meets  the  specifications  for  two  values  of  t,  i.e. ,  t  =  0  and 
T  “  T.  In  addition  y(t,  t)  for  O^^T  might  be  fixed  for  several  additional  values  of  t;  how¬ 
ever  for  the  purpose  of  this  exposition  two  will  suffice.  Obviously  if  y(t,  t)  could  be 
represented  by  relationships  similar  to  Equations  (67)  for  all  O^t^T,  the  specifications 
could  be  met  exactly.  The  number  of  values  of  t  for  which  y(t,T)  is  fixed,  for  instance  by  , 
Equations  (67),  will  also  be  the  number  of  terms  in  the  final  separable  function;  thus  the 
better  the  desired  approximation  the  higher  the  order  of  the  resultant  system. 
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Suppose  Figure  4  is  redrawn  as  in  Figure  5  in  which  the  t  axis  is  now  normal  to 
the  plane  of  the  paper.  Outside  the  region  O^t^T,  y(t,T)  is  arbitrary.  At  the  boundaries 
of  the  region  y(t,T)  satisfies  Equations  (67). 


i 

‘y(t,T) 

* 

y{tj)-y(t,o) 

* 

^.T) 

~1 

y(t,o) 

1 

y(i 

1 

APPROXIMATION  OF  y(t,  T) 
FIGURE  5 


If  the  only  constrai  ts  on  y(t,T)  are  at  t=0  and  t=T,  y(t,T)  can  be  approximated  in  several 
ways.  For  instance  let 

y(t,  t)  =  y(t,  0)  +|^y(t,  T)  -  y(t,  0)  u(t-T)  (69) 


where  u(t-T)  is  a  step  function  occurring  at  t=T.  y(t,T)  then  corresponds  to  the  full  line 
function  in  Figure  5.  If  Equation  (69)  is  Laplace  transformed  with  respect  to  t,  the 
resulting  transform  Y(t,s)  is  given  by: 


Y(t,s) 


-  y(t,o)  ,  y(t,T)-y(t,0)  ^  -sT 
s  s 


(70) 


e 


-sT 


Now  if - is  approximated  by  a  Fade  approximation  1,  for  instance. 


•sT 


^  2  (s-3/T) 


s  T 


i[^s^  +  4/T®+  6/T^ 


(71) 


then 


Y(t,  s)>  Y*(t.  s)  =  +  [  y(t,  T)  -  y(t,  0)]  - 


2 

T 


(s  -  3/T) 


s(s^  +  4/T®  +  6/T^) 


y  (72) 


The  inverse  transform  of  Y*(t,s)  will  be  denoted  y*(t,T).  Other  Fade  approximations 
could  be  used  in  Equation  (70).  Obviously  a  better  approximation  would  be  obtained  if  a 
higher  order  polynomial  were  used.  This  particular  approximation  was  chosen  because 
it  is  relatively  simple  and  because  the  numerator  is  of  lesser  degree  than  the  denominator. 
This  latter  condition  forces  the  second  term  of  y*(t,  t)  to  be  zero  at  t  =  0,  and,  hence, 
y*(t,  0)  =  y(t,  0).  Taking  the  inverse  transform  of  Y*(t,  s, ),  the  approximation  y*(t,T)  is 
obtained  as: 


y*(t,T)  =  y(t 


,0) 

m 


-2t/T  .  ,V2~t^  ,, 

e  sin  (  Y" 


+  y(t,  T) 


1  - 


,  -2t/T  .  ,V2”t_^ 

3e  sm  (  -Ip  +  ^) 


t  ^  T 


17 


1 


where  <jj  *  sin  ^  (1/3)  and  y(t,  0)  and  y(t,  T)  are  given  in  Equation  (67). 

Obviously  this  type  of  an  approximation  meets  the  differentiability  requirements  discussed 
in  the  previous  section.  As  a  very  rough  check  on  the  accuracy  of  the  approximation 
y*(t,T)  it  is  seen  that 

y*(t.  0)  =  y(t,  0) 

y*(t.  T)  =  0.  398y(t,  0)  +  0.  602y(t,  T)  (74) 

y*(t.  «)  =  y(t,  T) 

The  accuracy  of  this  type  of  approximation  can  be  improved  by:  (1)  approximating  y(t,T) 
by  several  steps  in  the  t  direction  (see  Figure  5),  i.e.,  let 

N 

y(t.T)  =  y(t.  0)  [y(t.  I  T)  -y(t>  (^))tJ  u(t-  g  T)  (75) 

—  sT 
N 

where  N  is  an  integer;  and  (2)  using  a  higher  order  Fade  approximation  for  e 

Another  method  of  approximation  which  might  be  used  is  to  approximate  a  higher 
order  derivative  of  y(t.T)  as  a  series  of  step  functions  in  the  t  direction  and  then  inte¬ 
grating  this  fvmction  an  appropriate  number  of  times  to  regain  y(t,  t).  As  an  example, 
the  derivative  of  the  dotted  line  representation  of  y(t,T)  in  Figure  5  is 


y(t,T) 


(76) 


If  this  equation  is  transformed,  the  Fade  approximation  in  Equation  (7)  is  used,  and  an 
inverse  transform  is  then  taken;  the  approximation  of  y^^'  ^^(t,T)  which  is  denoted 
y*  (t,T)  IS  given  by 


(77) 


where  again  =  sin  ^(-1).  Then  since 

y*(t,T)  =  y(t,0)  4-^  y*^°^\t,0)  dd 


(78) 


the  final  approximation  is: 


y*(t,T)  =  y(t,  0)“^'^^^  ||sin‘^+cos''^  j  +y(t,T)|l-e  sin^ 


+  cos 


'Cll... 


(79) 


y*(t„T)  =  0 


;  t<  T 
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As  a  comparison  with  the  approximation  of  Equation  (73)«  it  is  seen  that 


y*(t,  0)  =  y(t,  0) 

y*(t.  T)  =  0. 116y(t,  0)  +  0.  884  y(t,  T)  (80) 

y*(t,  «)  =  y(t.  T) 


This  procedure  then  produces  a  more  accurate  approximation  to  y(t. 


procedure.  Again  a  better  approximation  can  be  found  if  ^  | 

dr 

more  terms  as  in  the  case  for  y(t,  t)  in  the  first  procedure  (see  Equation  75). 


[y(t 


yt^.T 

,T)  J 


)  than  does  the  first 
is  approximated  by 


These  approximating  procedures,  while  leaving  something  to  be  desired  from  the 
standpoint  of  mathematical  sophistication,  are  straightforward  and  easily  applied  to  not 
only  weighting  functions  and  step  responses  but  also  polynomial  responses  which  are  also 
separable  functions .  The  advantage  of  this  method  is  that  the  functions  obtained  always 
have  enough  differentiability  to  apply  the  methods  developed  in  the  preceding  section. 

CONCLUSIONS 


In  this  paper  the  approximation  portion  of  a  synthesis  technique  for  linear,  non¬ 
stationary  feedback  systems  has  been  presented.  A  method  for  determining  a  weighting 
function  which  transforms  a  polynomial  input  into  an  output  which  can  be  described  by 
a  separable  function  of  two  variables  is  presented.  In  addition,  a  method  is  described 
by  means  of  which  a  function  of  two  variables  can  be  approximated  by  a  separable 
function  of  the  two  variables. 

In  a  companion  paper  18  the  synthesis  technique  is  completed  in  that  it  is  shown 
how  a  weighting  fvuiction  can  be  synthesized  as  a  feedback  system. 
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